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Abstract 



We study chaotic inflation driven by a real, massive, homogeneous minimally coupled 

scalar field in a flat Robertson-Walker spacetime. The semiclassical limit for gravity is 

_^ considered, whereas the scalar field is treated quantum mechanically by the technique of 

flj . invariants in order to also investigate the dynamics of the system for non-classical states of 

^^ ' the latter. An inflationary stage is found to be possible for a large set of initial quantum 

states, obviously including the coherent ones. States associated with a vanishing mean value 

of the field (such as the vacuum and the thermal) can also lead to inflation, however for 

p<J , such states we cannot make a definitive prediction due to the importance of higher order 

K^ ' corrections during inflation. The results for the above coupled system are described and 

0^ , their corrections evaluated perturbatively. 

O 

^ ! 1 Introduction 

t^ . 

Q>^ , A cosmological exponential expansion was proposed in the late 70's and early 80's in order 

"^ ' to solve the problems of the big bang model |^, g]. Many proposals have been made for 

^h| the realization of such an expansion, called inflation. The chaotic model proposed by Linde 

I ' 1^ was devoted to searching for conditions of a scalar field for which inflation may occur. 

f^£\] No restrictions, such as an initial state of thermal equilibrum or a minimum of an effective 

potential as in previous models ||, Q, were assumed. Indeed the strength of chaotic models 

is their relative insensitivity to initial conditions l5l . In this sense they are more generic than 

j^ ' other models, such as new inflation H| , which are more dependent on initial conditions and/or 

JH I require fine tuning of the parameters which enter in the potential of the scalar field. However, 

when chaotic models were first presented, it seemed that an implausible initial homogeneity 

was required in order for inflation to start. Therefore much effort was dedicated to showing, 

by numerical and/or analytical estimates, that anisotropics, small and large inhomogeneities 

will in general not suppress the onset of inflation (see the review [pj for a detailed discussion 

of this topic). 

In this work, limiting ourselves to the homogeneous case, we examine the idea of chaotic 
inflation from a semiclassical point of view, in order to enlarge the conditions which allow for 
inflation to arise from quantum mechanical initial states. Indeed, most of the chaotic models 
are also based on a classical behaviour for the scalar fleld. The " classicality" of the scalar 
fleld is a rather questionable assumption: let us note that the realization of an inflationary 
domain, usually as large as several horizon radii ||6| , implies that one is considering a region 
of space-time characterized by a typical size much smaller than the Compton wavelength 
(1/m) of the scalar field (to ^ i/ is usually considered in order to obtain sufficiently small 
density perturbations during inflation [@, Q|). On considering gravitational interactions it is 
difficult to think of the scalar field as being in a classical state when gravity effects it on 
distances smaller than its Compton wavelength. 

For the above reasons, it is of interest to study the evolution of the system for non 
classical states of the scalar field also. Here we present a quantum mechanical treatment of 



the homogeneous scalar field which drives a chaotic inflation model with a quadratic potential 
associated with its mass, V{(t)) = m? (/'^/2. This is greatly facilitated by the possibility of 
obtaining a Fock space of exact solutions to the Schrodinger equation for the matter field by 
solving the related harmonic oscillator problem with time dependent parameters, using the 
technique of invariants pO| . Our intention is to investigate whether quantum states of the 
scalar field exist which lead to a sufficient inflationary phase. 

In a certain sense we are analyzing " phenomenologically" the inflationary stage which is 
originated by a variety of possible quantum states. Therefore our conclusions could be of 
interest also for quantum cosmology, whose goal is the initial distribution of the scalar field 
after the quantum era lO . A classical initial distribution was considered in most of the cases 
treated in literature, except for some rare examples where quantum uncertainty was used to 
suggest a spread-out distribution for the classical value Q . 

We shall find an inflationary stage, sufficient to solve the problems of standard cosmology, 
also for states different from coherent ones ||lj, |lj] . In fact there is also the possibility that 
states such as the vacuum or a thermal ensemble (which have a zero expectation value for 
the scalar field) lead to an infiationary stage for which the time derivative of the Hubble 
parameter H is constant and negative. Indeed the time behaviour for _ff, which satisfies 
an interesting effective equation, is the same as previously found for the classical coupled 
system during inflation [|5|. However we are not able to study the evolution of the above 
non-classical states for arbitrary times because of quantum corrections (due to fluctuations) 
to the usual equations for matter. 

The above mentioned fluctuations arise naturally from the Born-Oppcnhcimcr (BO) re- 
duction pq , p^ , which we use to study the evolution of the matter-gravity system. In some 
cases these quantum corrections are not small and force the scalar field out of the initial 
quantum state. It is not clear whether this effect also alters the above mentioned qualitative 
behaviour of the scale factor during inflation. 

The outline of the paper is as follows: we set up our formalism in section 2 which is divided 
into four paragraphs, of which the first two treat the Wheeler-De Witt (WDW) equation [^ 
for the Robertson- Walker (RW) minisuperspace considered and its BO reduction. In the third 
paragraph we describe the method of invariants for time dependent harmonic oscillators and 
in the last we estimate the corrections to the usual matter and gravity equations of motions 
while relegating some useful formulae for their computation to the appendix. It is clear that 
in doing this we are taking seriously the WDW equation, which represents the full quantum 
version of this homogeneous cosmological toy model. The inclusion of all the Fourier modes 
of the scalar field would certainly represent a more reliable version of what is being studied 
here. However, we think that consideration of an exact quantum problem (although a finite 
dimensional one rather than an infinite), from which the semiclassical approximation can 
be obtained and verified is an important point worth studying. Section 3 is devoted to the 
study of the dynamical system. In the fourth section we present the numerical analysis and 
summarize our results in the Conclusions. 

2 Basic equations. 

2.1 The classical system 

We shall study one of the simplest cases of chaotic infiation, which is based on the theory of 
a real free massive scalar field (j) minimally coupled to gravity. The action is: 



S = d xC —Id X\J —g 



^-i..-5.0a.0-i^V 



(1) 



where L is the lagrangian density and m is the inverse Compton wavelength of the field 0. 
Further we consider a RW line element with flat spatial section: 

ds^ = g^^dx^'dx'' = -N^{t)dt^ + a^{t)dx^ (2) 

where N{t) is the lapse function and a{t) is the scale factor. On requiring homogeneity for 
the scalar field, we obtain through an Arnowitt-Deser-Misner decomposition, the following 



form of the action (|y) : 



S = I dtdxNa^ 



IGttG 



'(3) 



R 



K,jK'^ - K^ 




m^4)'^ 



(3) 



the extrinsic curvature tensor K^j and the extrinsic curvature scalar K for the metric hne 



where we have neglected boundary terms. The three- curvature scalar ^^'R, the square of 
the extrinsic curvature tenso 
element (0) are respectively: 

(3)i? = 0,if,,X^-' = ^^,X = -^-, (4) 



' '' ~N^a^' Na' 

where a dot denotes the derivative with respect to the time t. By using the relations 
the expression for the action (m we obtain the following lagrangian: 



L^ NV 






(5) 



where V is the volume of three space, M^j = 3/47rG and N plays the role of a Lagrange 
multiplier and is not a dynamical variable. The Hamiltonian is given by 



H„ 



NV 



2MI a 



2a3 



NV 



2MI a 



H 



M 



(6) 



where iTa = —M^iad/N and tt^ — a^(j)/N are the momenta conjugate to a and (f> respectively. 
The presence of A^ as a factor in (0) reveals that the total Hamiltonian has to be zero: this 
equation is usually called Hamiltonian constraint 



dHtot 
dN 







(7) 



and describes the dynamics, while the so called momentum constraints are trivially zero in 
this minisuperspace model. In order to quantize the system we rescale a — > aV^'^ absorbing 
the V factor into the Hamiltonian. Henceforth we shall work with an a having the canonical 
dimension of length. 

2.2 The quantum system and its BO reduction 

We start from a fully quantized gravity-matter system and subsequently return to a system 
in which gravity is semi-classical pjj{ and coupled to quantized matter, as is illustrated in 
||l|, |l^ . With a coherent state [^ |14[ for matter it is also possible to return to the completely 
classical case. 

On quantizing the Hamiltonian constraint one obtains the WDW equation |18| : 



Htot^(a,( 



92 1 , \ 



2A/p2 da"^ a 



(8) 



where ^ is a function of a and (j) and describes both gravity and matter. Subsequently one 
performs a BO decomposition of '5 as: 



*(a,0) = a7/;(a)x(a, 0) 



(9) 



where x{o-, 4') is not further separable. Let us add a brief comment on the quantization of the 
Hamiltonian constraint and on the assumption (y): in order to obtain eq. (jg) from eq. (0) a 

particular choice of factor ordering has been made i-;^ ^ ^^aa) q- The ordering ambiguity 
in this context has been previously noted |11[ . We observe that our choice is convenient since 
it is associated with eq. (g), which agrees with ^ = for a < as expected. All orderings 
will coincide in the classical limit and differ by 0{K). We have verified that such differences 



are negligible during the inflationary regime and in the semiclassical limit for gravity, which 
we consider. 

An effective equation of motion for ip may now be obtained p7[ by first substituting the 
above decomposition into eq. (g), contracting with x* and dividing by {x\x)- One then has 



2M2 



D^ 



(xI^mIx) 



(xlx) 



^ = - 



(xlD^lx) 



2M2 (xlx) 



V' 



2M2 (xlx) l^^'^a 



1 



where we have introduced covariant derivatives: 

iA; D 



da 



with: 



and a scalar product: 



A 



: ix\£x) 

^ (xlx) 



lx)(xh d 
(xlx) / da 



aa 



da 



Ix) W- 



(xlx) = / d0X*(a,0)x(a,' 



(10) 



(11) 



(12) 



(13) 



where the integral is over the different matter modes. 

An equation for x niay be obtained |lj by multiplying ( [lO|) by x and subtracting it from 
i): 



^^{ 



H, 



M 



(^M))x + T7r(D^)Dx 



A£ 



2^p/ 



^ 



da^ 



pi 
da^ 



V^ P' - (D')) X 



''da 



2Mp, 

da ^da' 



X- 



(14) 



We note that the r.h.s. of eqs. ( pi]| ) and (|lj) are related to fluctuations, that is they consist 
of an operator acting on a state minus its expectation value with respect to that state. 
Henceforth we shall refer to them as RHS fluctuations. 

It is convenient to write the above equations in terms of normal rather than covariant 
{D) derivatives and take into account the contribution from the connections by multiplying 
the gravity and matter wave functions by a phase. This redefinition leaves invariant the 
original wave-function '5 in (S) since the two phases are opposite in sign. With 



^ 



_ +i ("'Ada 



^ 



eq. ( |lO| ) becomes 



2M2 da^ 



a(H 



M) 



V^ 



'^K 



m^ 



(15) 



(16) 



Further we observe that in the semiclassical (WKB) limit for -if) one has: 



where Seff is solution to the Hamilton- Jacobi equation 



(17) 



1 



^Ki 



dS, 



eff 



do 



a{HM)=0, 



(18) 



and W = {dSeff/da) 2 . We analogously define x by 



X = e J 


A.a'^ 


so that eq. (14) can be written in the form 




a^{HM (HM))x+^j,^[gj) 


d , n^ 


9a ^ 2Mp2 


On introducing a time derivative 




^ d n 

in = —7. zr 


dSeff d 



^ {& - (D^)) x- 



and defining 



"dt ' M^i ada da 



Xs = e-i/*<-f^">''*'-''/°''°'-^;^ = e"^/*<-^^^>''*'x 



eq. (O) may be further rewritten as: 



H 



M 



ih- 



dt 



M^ia 



^-■k J*{HM)dt'-iJ'' da'A 
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(19) 



(20) 



(21) 



(22) 



X- (23) 



which is defined where ip has support in the scmiclassical limit (eq. (|l7|)) and henceforth we 
shall indicate by a the classical scale factor a{t). We note that on neglecting the rhs of eq. 
(|23| ) one has a time dependent Schrodinger equation for Xs- It is important to remember 
that with this approximation one does not take in account the transitions of the quantum 
state of the Schrodinger equation caused by the RHS. We shall discuss these transitions in 
sections 2.4 and 4. 

Let us stress that the peculiarities of this coupled system originated from the BO re- 
duction of the WDW equation. The Schrodinger equation, involving time defined in (pi|), 
describes the quantum evolution of matter, which in turn leads to the back-reaction appear- 
ing in the semiclassical evolution of the scale factor in eq. (p^. 

2.3 The quantum matter 

The quantum evolution of matter, as described by eq. (ESh, is given by a Hamiltonian: 



Hm — 



'2a3 9(^2 



(24) 



which is a harmonic oscillator with a time dependent mass a^{t) and a constant frequency m 
and has been previously studied (see for instance jlOl M ) without resorting to the adiabatic 
approximation pi] ]. 

However, since we are interested in a problem for which the evolution of a{t) could have 
an exponential behaviour, it is much better to employ a technique whereby a time-dependent 
problem can be solved exactly (we work in the Schrodinger representation). This technique 
is based on an auxiliary invariant operator (historically called adiabatic invariant |QQ|) which 
satisfies 



z?i— + [J,Hm]=0. 
at 



(25) 



The quadratic hermitean invariant originally introduced in pO| has real time-independent 
eigenvalues and is given by: 



+ (pTT^ - a^pc, 



(26) 



where p satisfies the foUowing equation: 



p + 3— p + rn p 



a' ' a^p'^ 



(27) 



This invariant / is important and useful since a general solution of the time dependent 
Schrodinger equation with the Hamiltonian ([24|) can be expressed as: 



x,o. = E 



c„e 



ifnit) 



n,t) 



(28) 



where the c„ are time-independent coefficients, \ n,t) are the eigenstates of the invariant 
(|2q ) with real eigenvalues A„ and the phases 9?„(i) are given by 



^n{t)^l I {n,t'\indt'-HM{t')\n,t')dt' 



(29) 



*0 



The quadratic hermitean invariant (Eq) can be decomposed in basic linear non-hermitean 
invariants lEOl: 



hit) = e»eW6(t) = ^ 



m 






^+1 [pn^-a" 



-ie 
/2h 



with 



0W = 



f-.( 



pn^- a-' 



The quadratic invariant (Gq) can then be written as: 



Iit) = h{llit)h{t) 



fi(b''{t)b{t) 



(30) 



(31) 



(32) 



with [/f,,/^] = [^: ^^] = 1 which allows us to identify the pairs of creation and destruction 
operators, {lb, ll) and (&, V'). The first is used to define a Fock space describing the solutions 
to the Schrodinger equation, and the second one generates the space of the eigenstates of 
the invariant defined in (Eq). As already shown in |22] this difference is important because 
X, as defined in (g), belongs to the second space, whereas Xs, as defined in (|2^), belongs to 
the first, of course always on neglecting the RHS fluctuations. 

Let us emphasize that the h operators that factorize the quadratic invariant are quite 
different from the creation and distruction operators (which we may denote by d) that 
factorize the Hamiltonian (p4|): 



with: 



Hm = Timid\t)d{t) + - 



(33) 



(34) 



and [d,d)] = 1. The two sets of d and h operators are related by the following Bogoliubov 
transformation 



B* {t)h + A* {t)P 



h 

2 



6t 



p{a^m)2 + 



p{a m) 



p{a^m) 



p{a^m)'- 



a 
ip\ — 



ip 



m 
m 



(35) 



and similarly for d^ . 

Since the quantum mechanical problem related to this time-dependent harmonic oscillator 
can be solved exactly, as in eq. ( p8| ) , let us first consider the state which best represents the 
classical behaviour of the scalar field in order to reproduce the previous results 1 23 , 24[ |5[ ^ . 
As is known |14t] this state is a coherent state psl for /;,, which satisfies: 



It\ a,t)^ = a\a,t)^ (36) 

where a = \a\ e^^ is a time-independent complex c-number. Such a state satisfies the 
Schrodinger equation and can also be written as a normalized superposition of eigenstates 
of the quadratic invariant in (E6f): 






°° n -inQ 



n=0 



E 

n=0 



n,t] 



(37) 



On defining (O)^ = a{oi,t \0\ a,t) ^, the mean value of some physically relevant quantities 
is: 



7r0,c 



>0>a = V2h\a\ 



'^\a\pcos{e- (3) 



a^ p cos(e -13)-- sin(e - /?) 



(38) 
(39) 



((A<^)2)^^((^-0,)2)^^_p2 



(40) 



; (A0)2 )l ( (A,r^)2 )i = I yiT^V^ 



We note that the mean value of the field (38) satisfies the classical equation of motion: 



+ 3 — (pc + m (/)c = , 



(41) 



and a does not appear in the uncertainty relation {Mm , which however increases with time. 
The mean value of the matter Hamiltonian for the coherent state is: 



H, 



M , 



>0>'+"^V 



1 

Hci + Ho = a^ipci +Mo) 



l + ^{p'a^ + ^ + m^p'a' 



(42) 



where Hd is the classical Hamiltonian for the scalar field (f>c defined in eq. (Bq) , Hq is the zero- 
point Hamiltonian and fici, fJ-o are respectively the classical and zero point energy-densities. 
All these quantities are time-dependent. In order to consider the above semiclassical field as 
a fluid we further define, in analogy to the classical case, a pressure operator: 



1 



2„3I2 



p = — ^ -^ - rira-^(j)- 



2a3 



a" 



(43) 



Its expectation value in a coherent state is: 

^2 * 



yP)o.= 



1 
2^ 



0:C 2 3i2 I ' / -2 3 

^ - TO a-'cj)^ + 2 i ^ "■ 



a^ p^ 



m^ p^ a^ 



: Pel + Po 



and we have the usual fluid conservation law: 



(44) 



Ac/ + Ao + 3- {p-cl + P0+ Pel + Po) = . 
a 



(45) 



2.4 Higher order corrections to semiclassical theory 

We are interested in the corrections given by the rhs of eq. (p^ for gravity and by the rhs 



26| , for which (x|x) = 1- For this 
) can be constructed by using the 



of eq. (|1J) for x in a coherent state or an n-particle state 

purpose we note that the x wave function introduced in _ 

b, b^ operators. Let us then indicate by | a)^ the coherent state associated with the operator 

b while the | n)^ are the eigenstates of the associated number operator. Further, we shah 

derive expressions vahd during the inflationary regime (a >> 1), so that sonic simplification 

will occur. 

Let us first consider the corrections to the gravity equation. With the introduction of the 
semiclassical time in eq. (pi|), we can write: 

On defining the Hubble parameter as iJ = ^, the semiclassical version of eq. dlO) leads to: 

^' = ^ix\HM\x) + ^Xmx)-^{HM) + ^;,^^P.. (47) 

With the help of formulae in the appendix one obtains: 

Pn = ||| n)f - \,{n ||| n)f = 2m'\A\'\B\^{n^ + n + I) (48) 



d . , „ , , .d 



'dV "" '"' 'dV 



Pc.^\^Jc^)f-\t{a\^\a)f ^ e-m(l + 2|Ap) \a\' + 2mM^\Bmi + 2\a 



2 



2m {a^AB* + a* ^BA*) (m (l + 2\A\^) - 9) . (49) 



We note that for the coherent case the terms in a of order higher than quadratic cancel. For 
completeness we also write the average matter Hamiltonian for an | n)i^ state: 

b( n \Hm\ n), = hm{\A\^ + ^){2n + 1) = aVo(2n + 1) (50) 

and on solving eq. (p7[) for H^ one obtains 



•l + ,/l ■ '''^^ 



Let us now estimate the corrections, for a 3> 1, to the results obtained when the rhs is 
neglected for the coherent, the | n)j, and the vacuum states respectively. One first observes 
that the time dependence of the rhs in this limit is the same as that of the Hamiltonian, so 
that their ratio becomes time independent. For the coherent case in the |a| ^ oo limit we 

have: 

h^Pc, 2 + cos 2/? 

2 ~ .I„I2 (52) 



b 



'a\HM\a)l 4|a 



which goes to zero, so that one may expect no correction in the classical limit. For the \ n)^^ 
state case, in the n ^ oo limit, we have 



h^Pn _ , 1 

2 



2 



b{n \Hm\ n)j 



(53) 



which leads to an 11% correction. Finally for the vacuum case (a = n = 0) we obtain: 

,2-2 (54) 



n^Po 



b 



{o\HAi\on 



which gives a correction (^ — 1) of order 36%. 

We now consider the matter equation. It is convenient to introduce 

I '/'X2,xi> = I X2) - (Xi 1X2)1 Xi)- 



(55) 



As the corrections are an operator we shall examine their matrix elements with respect to 



the 



states and the coherent states (an overcomplete non-orthogonal set). The general 



will be 


(^x.,xl^^lx)+2(xl^lx)(0x..xl^lx)J 


2M>a2 



(56) 



In the large a limit, during inflation, the first two terms of eq. ( p6|) will be subdominant 
with respect the last two by a factor of order a^, hence we may neglect them and denote the 
last two by Qx2,x- 

Using the results given in the appendix and in eqs. (|5y) and (|5^) one has for the ratio of 
the the matrix elements between the | n)j, and \ rn)f^ states and the average Hamiltonian : 



^va.7i 



1 



-4,\Jn{n~ l)6m,n-2 + \/n{n - l){n - 2){n ~ 3)6„ 



4 + 



b{n\HM\n), 4e(2n+l)2 

v/(« + l){n + 2){n + 3)(n + 4)S,n,n+4 + W{n + l)(n + 2)S„^^r.+2 



(57) 



which are of order unity for the vacuum {n = 0) state and one order smaller, but finite, 
for large n. The form obtained reveals that the allowed transitions will eventually lead the 
system to spread over a superposition of states 

In the same manner we obtain expressions for the coherent states in the large a limit. 
Only the last two terms of eq. {b&j contribute giving 



b{f3\a)bn^m^\A\' 
2M2a3iJ2 



iP* - a*) {a* +p* + 4a)(4 + p*^) 
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3/3* 



4aa* 



(58) 



In order to see how large such corrections are in the |a| — > 00 limit, when the system is 
characterized by a behaviour close to the classical, we set 7 = /3 — a and obtain 



Qp,a 



> \Hm\ a)b 



4£.b{a\HM\a)l 



8a* + 16a + 47* + 3(c 



+ 7 a + 



7*^ + 8aa*^ + 8|ap7* + 4a7*^ + 8a^ + 14a*a^ 



|q|>>i |a| 



(59) 



We note that the above quantity is bounded in 7 so that in the |a| — > cx) limit every transition 
away from a coherent state is suppressed and thus, as expected, one recovers the classical 
behaviour of the system. Further the behaviour 0{-^) is also valid for any a. 
Let us finally exhibit the ratio for the vacuum case (a == 0): 



0/3,1 



\0\Hm\Q 



46(0|i?Af|0)^ 



(/?*■* + 4/3*2) 



4? 



(/3* 



.4/?*2)e-l^l /2. (60) 



Of course if one multiplies the above by b( n |/3)b and integrates J d/i(/3), where the measure 
(d/^(/3) = d.ne0dimfi -^ jg g^^j^ ^^^^ ^^g completeness relation is / d^(/3)| /3)j,6(/3 | = 1, we 
re-obtain the matrix elements for n = in eq. (m) as expected. 

Since we are also interested in the time evolution for the system initially in the vacuum 
state we shall also estimate in section 4 the transition rates (to the | 2)^ and 14)^^ states) 
which are generated by the correction terms on the rhs of the matter equation. 



3 The analysis of the dynamical systems 

In this section, in preparation for the numerical analysis, we first derive an effective equation 
describing the evolution of semiclassical gravity after "averaging" over the quantum matter 
degree of freedom. In our model, on neglecting the RHS fluctuations, this equation will 
not depend on the matter quantum state. Subsequently we shall study the dynamics of 
the matter quantum state in order to understand and possibly evaluate the effect of the 
approximations made. In particular we shall consider in detail the evolution for a coherent 
state, a vacuum state and a thermal state. 

3.1 Effective equation for gravitation. 

Let us begin by considering matter in a generic quantum state \ x) s satisfying the Schrodinger 
equation with Hamiltonian H]\i . The backreaction on gravity is described by the equation 

H' = -^^six \Hm\ x)s = -^six lAl X>. = ^-M (61) 

where fi — Hujc? is the energy density operator and /i its average. On differentiating 
equation ( |6l|) with respect to time one obtains: 

^^\^^^ ^^^^ ':l^' ^^^ = -4^G,(x IpI x), ^ -4^Gp (62) 

where p is the pressure operator defined in (B3) and p its average. One notes that ij = 
— 47rG(p + /i) so that a fluid-like equation is satisfied: 

/i + 3iJ(p + /i) = (63) 



We now differentiate eq. (|6^) with respect to time and use the relations ( |61| ) and (|62 
in order to eliminate the average, with respect to | x) si °f ^^ Hamiltonian and of its time 
derivative. One then obtains 

ii^ + 6i/i/=^m^(x|{0,Mlx). (64) 

where the term on the RHS comes from the commutator [p, iJ^i/]. Further we perform yet 
another time derivative of d64) obtaining 

1 97rG 4?7rG 

H+6H' + 6HH = -—^m'Hsix \{<l>,^.p}\ x), + -^^'^ix | [i?M, {^,7r4]| x), (65) 

One easily computes the commutator in the last term [Hm,{4'tT^cI>}] — —4:iha^p, then, on 
using the relations (p2|) and (p4), one finally obtains an equation for H only: 

H+6H^+4:H + 9HH + 18H^{H+l):^Q, (66) 

o 

where we have scaled the time by m^^ and henceforth we shall use this new scaled time. 

The above equation is valid for every quantum state, in particular for a coherent state, 
and, on taking the classical limit for matter, also for the classical field. Thus the classical 
model of chaotic inflation can be described by this equation. Indeed, it is known that the 
classical system is governed by the system of equations 



^2^ SttG 



(^) 



4 + 3*.\/M(*l±i£!)+,. = (67) 

which gives a non linear second order differential equation in H but, as is easy to check, 
also implies eq. (p6h. Hence the classical system has a dynamics which lives in a two- 
dimensional surface in the (H, H, H) phase space. This surface is an invariant manifold 
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Figure 1: In gray we have the invariant manifold corresponding to the classical dynamics (|67|). 
Also exhibited are a trajectory related to a highly non-classical state with initial conditions 
{H = A, H = 0.4, H = 0) and a line, given by {H = —1/3, H = 0) and not intersecting the 
invariant surface, which is approached by the system in the inflationary phase both for classical 
and quantum matter; after reaching it the trajectory will lie close to the invariant surface. Let 
us note that in the figures (1 and 2) we use a prime to denote a time derivative. 



under the dynamics of eq. (|66| ) which was obtained for quantum matter. It is the initial 
conditions that distinguish the two cases: in the classical case there is a precise relation 
between H{to),H{to) and H{to) which is not satisfied in the quantum case as can easily be 
seen, for example, for the vacuum or a coherent state. The trajectories for a generic quantum 
state do not lie on this surface, but they approach it asymptotically during evolution, however 
they will never exactly lie on it (see Fig. 1). 

In order to understand how the inflationary regime can be reached one may study the 
set of initial conditions in the three dimensional {H,H,H) phase space which lead to it. Let 
us note that the only fixed point in the finite region for the vector field associated with eq. 
(pq ) is given by {H, H, H) = (0, 0, 0). This point is asymptotically stable and the divergence 
of the vector field is negative everywhere for positive iJ, so that every trajectory will fall 
towards it. 

On the basis of the numerous analyses already performed on the classical system [n5l 



24 



lb\ one expects the value of H to change during the eventual inflationary phase, 
decreasing (on starting from initial conditions on the two-dimensional invariant manifold 
shown on Fig. 1) from values generally much greater than one in the units we use (Planck 
mass). This implies, on comparing the magnitude of the different terms of the eq. (p9), that 
one can expect, from the last two terms, unless H is small. 



TJwO 



H 



(68) 



Let us note that this is the same behaviour as the classical case, and is quite different from a 
true de Sitter evolution for which H = const. We remark that this is the same behaviour as is 
also found by the authors of Q, who describe the curve corresponding to (g8|) by 0c ~ const. 
Therefore in the three dimensional phase space, the trajectories starting from different points 
( which are related to the initial conditions), for both the classical and quantum matter cases 
will, if leading to an inflationary regime and after an eventually transient phase, converge to 
a region around the curve characterized by (pSl) and a decreasing H. This is what is shown 
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and commented in Fig. 1. Since one is also interested in the pressure/energy-density ratio, 
let us write the expression to be used later in the numerical simulations: 

p 2H 

remembering that it is valid for any quantum state satisfying the Schrodinger equation. 

The above approach for gravitation is general and does not require any direct solution of 
the Schrodinger equation for matter. However it is important to know the matter state, and 
for this reason we introduced the technique of the invariants. On using them it is possible 
to construct coherent states to study the classical limit and estimate if a quantum state, 
satisfying the Schodinger equation, defined by the Hamiltonian Hm will be modified by RHS 
fluctuations, possibly before the end of inflation. We shall attempt perturbative estimates 
in section 4. 

3.2 Matter states. 

For the coherent state one has a four dimensional dynamical system, which is given by eq. 
(1^) for p and: 

(70) 
(71) 

We are interested in evolution from an initial state which is well classified in terms of the 
Hamiltonian at the initial time (t — to), that is a coherent state constructed from the 
particle annihilation and creation operators which factorize the Hamiltonian. This choice 
corresponds to requiring the invariant and the Hamiltonian coincide at t — to and the 
Bogoliubov coefficients will be A{to) = and B{to) = 1, leading to the conditions: 



a) 


= ^^(/ic/ +A*o) 




© = -r^ 

a-^/9^ 



p\t 



1 







TOa^(to) 
Kh) = (72) 

Subsequently we shall consider a larger set of initial states not satisfying the above. It is 
remarkable that the above requirement relates the fluctuations of the scalar field given by 
eq. (|4^) to the scale factor of the universe. Indeed at the quantum level the scale factor 
a appears inextricably in the equations (already in the energy density, see for example eq. 
(E3)), this is the basic reason for which the semiclassical system is three dimensional. From 
eq. ([T^, corresponding to a possible choice of initial conditions, it is clear that the smaller is 
the scale factor the larger is p, which would actually be infinite if a were zero as is classically 
allowed. 

Let us now examine the vacuum state. This state is interesting even if it is obvious that 
it will be the first to be changed by the transitions generated by the RHS fiuctuations of 



the matter equation ( |23[ ). The vacuum case corresponds to a = {p,ci = in eq. (|70|)), so 



that O decouples and we may eliminate a in eq. (|27| ) obtaining an equation for p and H . In 



analogy with eq. (62), on scaling p by vGh, one also obtains: 



•2 _ J_jj2 

p + 3Hp + ^ + 2p = (73) 

P 

H + 3H^ - 27rp2 = . (74) 

In this case we study the evolution of an initial state annihilated by the invariant destruction 
operator which is not necessarily the Hamiltonian one, unless the conditions given in (|72) 
are imposed, and parametrize it in terms of the initial conditions (p(to), p{to) and H{to)) or 
{H{to), H{to) and H{to)). We observe that if the scalar field is in the n eigenstate of the 
invariant , the dynamical system is given by eq. ( [73| ) with the term H^ divided by (2n + 1) 
and by eq. (|74| ) with the term p^ multiplied by {2n + 1). 
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We may now consider a thermal (not pure) initial state which is also characterized by 
zero mean values of the matter field and of its derivative. Such a state could arise within 
the context of our BO approach on introducing a complex, rather than real, time, thus 
generalizing our previous considerations. In order to describe such an initial state we consider 
a density matrix constructed using as basis a given number of quanta of the Hamiltonian (at 
t = to)- If at time to the conditions given in ( |7^ ) are satisfied, then the time evolution of the 
basis will be given exactly by the eigenstates of the quadratic invariant so that the density 
matrix will be 

p^(l-e-^")^e-'5"'=|fc,t)(fc,t| (75) 

fc 

where /3 = l/ksT and of course during evolution it will no longer be a thermal state of the 
Hamiltonian. The operator p satisfies the Liouville equation so that also for this state one 
easily derives the effective equation (|66| ) . In fact using the known evolution equation for the 
p operator we can compute, on taking traces, the mean value of all interesting operators. 
In particular one has for the average matter Hamiltonian density: 

r. (,*!)., .(1 + -^) (76, 

and for the matter fluctuations 

Tr{P^')-\p'{l + ^^)- (77) 

These differ from the vacuum case by an extra factor which increases with the temperature, 
as expected. Thus, even for not so large values of p we could still have large fluctuations of 
the matter field if the temperature is high. 

4 Numerical Simulations. 

In this section we exhibit the results of the numerical analyses of the previously described 
dynamical situations. 

Let us first consider the basic effective equation (p6|), which is valid for every quantum 
system within the framework of the approximations made. We find that, starting from initial 
conditions very far from the "classical" surface, the inflationary regime coincides with the 
region described by eq. (|68| ) and the trajectory will lie close to the two-dimensional invariant 
surface previously mentioned. This is illustrated in Fig. 1. 

In order to have a better idea of the different possible trajectories arising from differing 
initial conditions we show in Fig. 2 a projection of the trajectories on the H-H plane. From 
the picture it is evident that trajectories will lie, to within a very good approximation, for 
some time on the curve given by (pq). Therefore it is evident that we have a large set of 
initial conditions which will lead to an inflationary stage. One also sees that the number 
of e-folds reached at the end of the inflationary stage depends on the value of H at the 
moment when each trajectory approaches this curve. In fact it can be seen that during the 
inflationary stage starting at i = t, one has 

^«exp(^i/(iOi + i7(t,)y) (78) 

and thus the total number of e-folds is 

N= f' H{t)dt^-^0^ (79) 

where tf denotes the time at which inflation stops and we have considered H{tf) <C H{ti). 
In order to compare our results with the already known behaviour of classical models 
|15| , p^ , p4| , pSJI we shall start from a system described by a coherent state which will be 
squeezed during time evolution. The results of the analysis are easy to understand from the 
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Figure 2: Trajectories with different initial conditions projected on the H-H plane. The bold 
lines reflect a superpositions of curves having quite different H for some time. It is evident that 
all the curves reach the common region described by eq. ([6^). 



structure of equation (pffl) on noting that H (or a) is completely governed by the matter 
Hamiltonian density which has two contributions, one from the vacuum and another which 
is exactly that expected from classical dynamics. 

Therefore if we are in a coherent state characterized by an average number of quanta 
|ap 3> 1, the vacuum contribution will be negligible and the resulting evolution will be close 
to the one given by the classical homogeneous matter field. For such a case we show in Fig. 
3 the e-fold number, the Hubble parameter, the mean field fluctuation, the pressure/energy- 
density ratio and the average field (a \4>\ a), which corresponds to the classical quantity (fic, 
resulting from the dynamics. We start from a condition, known to lead to inflation, with 
the average value of the field starting the slow rollover phase; more precisely, working in 
units of the Planck mass, we have taken for the inverse Compton wavelength m = 10^^, 
imposed initial condition (|2|) and |a| = 2.334 10^, /3 = 0, p{to) = 10"^, p(to) = so 
that (j){to) = 3.3, (pito) = 0. The resulting initial conditions for equation ( p6|) are H — 
6.7539, H = -6.283 10"^ and H = 6.878 lO^^^ (all these numerical values are for the 
rescaled dimensionless quantities and time scaled by m~^). After a short transient, a linear 
decrease of H and p until the inflation period ends (with an e-fold number of 70) is evident. 
The behaviour of the ratio between the pressure and the energy density remains, as expected, 
close to a value of —1 for most of the inflationary period and later starts to oscillate. 

From the point of view of the matter state, it is interesting to note that the role played 
in the classical case by the value of the field is now replaced by p, which is related to the 
quantum fiuctuations of the homogeneous field. During inflation p decreases linearly and 
after this "dissipation" it will start to oscillate above zero. We also note that the average 
value of the field behaves as the corresponding classical quantity, decreasing during inflation 
and oscillating after it (in the figure we use a logarithmic scale which covers only positive 
values) . 

Let us consider the other extreme case: that of the invariant vacuum (which eventually 
can also be the Hamiltonian vacuum at t = to). For this state the average field and time 
derivative of the field are zero. On neglecting the transitions induced by the RHS of the 
matter equation, it is clear that for high enough fluctuations (large values of p) the system 
could undergo an inflationary phase. In such a case the behaviour of a, H, p, p/ p is the 
same as that of the coherent case plotted in Fig. 3. However the presence of a non negligible 
term on RHS implies that the Schrodinger equation will be a rough approximation to the 
real matter equation and we are then unable to predict the evolution of the system. 

Since the vacuum is an extreme case we feel that it is interesting to try to further analyze 
such a situation. Let us then estimate perturbatively the time necessary for the corrections 
to produce a large transition rate such that it would be meaningless to consider the normal 
(Schrodinger) evolution for the vacuum. After such a time we shall no longer be able to 
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Figure 3: For matter in a coherent state we show the time behaviour of the e-fold number, the 
Hubble parameter H, the mean value p for the fluctuations of the field, the pressure/energy- 
density ratio and of the average value of the field (pc = {(f>)a- 
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Figure 4: For an inflationary region with different values for p (for H) which would lead, in 
absence of RHS correction, to an inflation with e-fold number of 70, 25, 15 , the corresponding 
(from left to right) time behaviour of the transition rate eq. (|8l| ) is shown. For the above the 
transition rate becomes significant and possibly modifies the inflationary regime after 1%, 4% 
and 6% respectively of the inflationary period in the absence of RHS fluctuation corrections. 

make any prediction and if the system was, or was about to enter, in an inflationary stage 
it may or may not thus continue. We rewrite the full matter equation (with tlie RHS) as 

(i?i^-i?M)|x). = -e^^i^xlx) (80) 



where (/? and F^ are defined by comparison with eqs. (03) and (^3|). The full higher order 
non-linear equation is clearly too complicated to solve. Hence we estimate the transitions to 
other states which are solutions of the unperturbed Schrodinger equation, as is customarily 
done in perturbation theory. This means we expand \x)s ~ Era'^'^(^)l ^) s ^^"^ evaluate 
the contributions of the perturbation to lowest order pTl , thus we consider the operator i^o 
which depends on the | ) state. One then obtains for the transition rate: 

-i?i^ = e™«(n|Fo|0) (81) 

Further we must choose the initial values for p, p and a. The following relations, obtained 
from eqs. (|4|) and ( |70|) are useful for this purpose: 



We first consider the case for which the state is initially in an inflationary regime so that 
the relations in ( p8|) are satisfied. On requiring a certain number of e-folds we find the value 
of H needed and hence the necessary p and p and a. For example values H > 1 will lead to 
p « 0.7H and p « —0.23 in the chosen units. On using these values we have computed the 
quantity in eq. (B^). In particular in Fig. 4 we show the real part of C4 for three different 
cases leading to an e-fold of 70, 25 and 15 if the corrections were not present (the right hand 
side of Eq. (Bw has comparable real and imaginary parts, for this reason the contribution 
to the real part of C4, plotted in Fig. 4, due to a possible rapidly oscillating phase of C4 
is excluded). We note that for smaller initial H (which would lead to less expansion) the 
transition rate becomes significant at a later time. If, on the other hand, we consider an 
initial vacuum satisfying the conditions of eq. d72) (as yet we are not in an inflationary 
regime), the transition rate is non- negligible from the beginning and as a consequence we 
simply cannot use the Schrodinger equation to evolve it (analogous considerations hold for 

C2). 
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Little can be added on the evolution from a thermal state, on neglecting the RHS, since 
it is similar the vacuum case. The basic difference is that the thermal factor in eq. (um will 
increase, depending on the temperature, the energy density. 

Let us end this section by adding a brief comment about reheating. It is known that an 
inflationary stage during which the scale factor is accelerated should be associated with a 
massive entropy production in order to solve the problems of the standard cosmology p8| . 
However, also for states of the matter field with a vanishing mean value persisting till the 
end of an inflationary phase, p would oscillate with the same period the classical field would 
have had, but its amplitude would decrease in a more complicated way (see, for instance, 
eq. (uSh for the vacuum case). The similarity of the behaviour of scalar field fluctuations 
for classical and vacuum states has already been noted [M in an exact de Sitter space-time 
for long wavelength modes. For these reasons in a model where a second field x is coupled 
through a 0^x^ term to the inflaton, the vanishing value of the scalar field does not seem to 
prevent the possibility of a reheating process after inflation. In any case we realize that the 
inclusion of all the other modes is necessary to tackle this program ^^ . 

5 Conclusions 

In this paper we have studied a massive homogeneous minimally coupled scalar field in a flat 
RW metric. After having quantized the system we recovered, through a BO decomposition, 
the semiclassical evolution for gravity and a Schrodinger equation for matter, always on 
neglecting RHS fluctuation corrections. Such a treatment is however an approximation which 
could break down for a matter quantum state because of the growth of the corrections. Since 
we do not know how to predict exactly the dynamical consequences of the RHS fluctuations, 
we evaluated them perturbatively and we studied in detail the properties of the dynamical 
system and the time scale at which the corrections are no longer negligible. 
In particular wc found: 

(1) The dynamical semiclassical gravity system for a generic quantum state of the scalar 
field is described by a non-linear third order differential equation for the Hubble parameter. 
We have seen that the classical system lives in an invariant surface of the three dimensional 
phase space and that an inflationary regime (with H ~ — i) is experienced for a generic 
quantum state of the scalar field. 

(2) The time scale at which the corrections are no longer negligile depends on the "clas- 
sicality" of the quantum state for the scalar field. In the classical limit (|ap -^ -l-oo) these 
corrections are suppressed while for a generic quantum state they are no longer negligible 
before inflation ends. For the vacuum state the corrections could become important, in the 
least favorable case (for example by choosing particular initial conditions), at the same time 
scale as that for which inflation starts. 

In any case we feel that from this analysis it emerges that inflation can arise from a quantum 
mechanical initial matter state, or at least a coherent one, and not just a classical state. 
Such a situation was envisaged in the earlier pioneering works [Q, pl[ . 

Let us conclude by again noting that the main problems encountered in considering 
inflation from a non-classical state of the scalar field are related to the difficulties in solving 
the full matter equation (|8^). One would like to at least evaluate the effects of the RHS 
fluctuations non-perturbatively and see how an inflationary regime evolves from a generic 
quantum matter state. We intend to investigate this and other points in more detail, as well 
as applications of the above method to other physical problems [p2i . 
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7 Appendix 

We here compute the time derivatives of the |n)f, and of the coherent states for the destruction 
operator b satisfying 



The coherent state | a)j, is defined by applying the operator D 



(83) 



y.b^ —a.^'b 



e 2 e e 



ab' ^ — Oi b 



D = e 

to the vacuum state | 0)^, for which b\ 0)f^ — Q. Therefore 



(84) 



dt 


.-,„,,,^..-. 


n=0 " 




^1 \ 


dt dt 


|0)fc 



with 



where, on recalKng the equations satisfied by Xi we have for the second term in (|8 



dt 



On using 



and 
we have 



»{0|^|O)j-imB/l-it^|0).. 



f 4|it,//„| + .eM^ 



[b\HM] = hm[-{l + 2|yl|2) fet _ 2ylS*&] 



|-fet» = j[e - m(l + 2|An]n6^" - 2^m5M[ "^'^ ^) fct"-2 + ^fn-i^] 



One finaUy obtains: 
d 



d 



(85) 



(86) 



(87) 



(89) 



(90) 



-\n), =. {zn[e-mi2\A\' + l)]+,{0\-\0),}\n), 

-imB*Ay/n{n-l)\ n - 2)^ - imBA* y^ (n + l){n + 2)\ n + 2)^ (91) 



and 
^^\a),^l-imBA*P' 

where 



ie - im (l + 2\A\'^ 



6t 



d_ 

'at' 



b^O 1^1 0)^-imAS*a2 



b(0\■|-^\0), = -^^o-^mi\A\^ + ^) 



(92) 
(93) 
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